Field-Driven Hysteresis of the d=3 Ising Spin Glass: Hard-Spin Mean-Field Theory 
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Hysteresis loops are obtained in the Ising spin-glass phase in d — 3, using frustration-conserving 
hard-spin mean-field theory. The system is driven by a time-dependent random magnetic field Hq 
that is conjugate to the spin-glass order Q, yielding a field-driven first-order phase transition through 
the spin-glass phase. The hysteresis loop area A of the Q — Hq curve scales with respect to the 
sweep rate h of magnetic field as A — Ao ~ h b . In the spin- glass and random-bond ferromagnetic 
phases, the sweep-rate scaling exponent 6 changes with temperature T, but appears not to change 
with antiferromagnetic bond concentration p. By contrast, in the pure ferromagnetic phase, b does 
not depend on T and has a sharply different value than in the two other phases. 
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Frustration and non-equilibrium effects induce compli- 
cated ordering behaviors that challenge the methods of 
statistical physics. Perhaps the most ubiquitous non- 
equilibrium effect, hysteresis is the current topic of in- 
tense fundamental and applied studies. The effects of 
a spatially uniform linearly driven magnetic field on the 
ferromagnetic phase of the n-component vector model in 
d = 3 |l[ and of the pure (no quenched randomness) Ising 
model in d = 2, 3, 4 0, y] have been studied. In these 
works, the scaling exponents of the hysteresis loop area 
with respect to the sweep rate are obtained. Spatially 
uniform pulsed, sinusoidally oscillating, and stochasti- 
cally varying magnetic fields on the ferromagnetic phase 
of the pure Ising model in d — 2, 3 and spatially uni- 
form square-wave magnetic fields on the ferromagnetic 
phase of the pure Ising model in d = 2 [8j have been 
studied. Two current detailed experimental studies of 
hysteresis are in 
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FIG. 2: (a) Equilibrium spin-glass order parameter Q ( ' as a 
function of temperature T — J -1 . The curves, from top to 
bottom, are for p = 0, 0.1, 0.2, 0.3, 0.5. The latter two curves 
overlap, (b) Equilibrium magnetization A/' ' as a function 
of concentration p. The curves, from top to bottom, are for 
T = 0.5, 1.0, 1.5,2.0,2.5,3.0. 



FIG. 1: Phase diagram from hard-spin mean-field theory for 
the d=3 Ising spin glass. All phase boundaries are second 
order. 

In the present study, hard-spin mean-field theory, de- 
veloped specifically to respect frustration [i| [l(| , is used 



to study the non-equilibrium behavior of the field-driven 
first-order phase transition that is implicit, but to-date 
unstudied, in spin-glass ordering. For the Ising spin- 
glass on a cubic lattice, the phase diagram is obtained 
and the temperature- and concentration-dependent or- 
dering of the spin-glass phase is microscopically deter- 
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FIG. 3: Zero-temperature spin-glass order parameter Q^- ' 
as a function of antiferromagnetic bond concentration p, ob- 
tained by averaging over 10 realizations, with the standard 
deviation being used as the error bar. The dashed line indi- 
cates the transition between the two phases, whose position 
is obtained from the phase diagram in Fig.l. 



mined. The random magnetic field that is conjugate to 
this microscopic order is then identified and used to in- 
duce a first-order transition and hysteresis loops. We find 
qualitatively and quantitatively contrasting scaling be- 
haviors in spin-glass, quenched random-bond ferromag- 
netic, and pure ferromagnetic phases of the system. 
The model is defined by the Hamiltonian 



(1) 



where Sj = ±1 at each site i of a cubic lattice and < ij > 
denotes summation over nearest- neighbor pairs. The 
bond strengths Jij are equal to — J with quenched prob- 
ability p and +J with probability 1 — p, respectively cor- 
responding to antiferromagnetic and ferromagnetic cou- 
pling. Hi(t) is a linearly swept quenched random mag- 
netic field, itself determined, as explained below, by the 
spin-glass local order of this system. 

For our calculations we use the hard-spin mean-field 

theory d El El El El El El El El El Hi HI m [22 , 

[23l | , a method which is nearly as simply implemented as 
the conventional mean-field theory but which conserves 
frustration by incorporating the effect of the full magni- 
tude of each spin. The self-consistent equation for local 
magnetizations m, in hard-spin mean-field theory is 
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where the sum {sj} is over all interacting neighbor con- 
figurations and the sum and the product over j are 
over all sites that are coupled to site % by interaction 
Jij. The single-site probability distribution P(rrij,Sj) 
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FIG. 4: The top-row figures are from a hysteresis loop in 
the ferromagnetic phase with quenched random antiferromag- 
netic bonds, T = 1.5, p = 0.15, h = 0.005. The middle-row 
figures are from a hysteresis loop in the spin-glass phase, 
T = 1.5, p = 0.4, h — 0.005. Left: calculated equilibrium 
local magnetizations m- ' in a cross-section of the three- 
dimensional system. A hysteresis loop is started from these 
systems. Middle: local magnetizations m;(t) at the first can- 
cellation point, M(t) = (top row) and Q(t) = (middle 
row), of the first hysteresis loop. Left: local magnetizations 
at the first reversal point, M(t) = -M(0) and Q(t) = -<3(0), 
which occurs when the first hysteresis half-loop is completed. 
The bottom cross-section shows the vanishing equilibrium lo- 
cal magnetizations everywhere in the paramagnetic phase, to 
be contrasted with the spin-glass cross-section immediately 
above it: the global magnetization M 1 - ' = in both cases. 



been used in time-dependent systems, in the study of 
field-cooled and zero-field cooled magnetizations in spin 
glasses. |17j| 

Equilibrium Phase Diagram - The equilibrium local mag- 
netizations mf^ are determined by simultaneously solv- 
ing N coupled Eqs.([2]) for all N sites i of the system, at 
zero external magnetic field, H = 0. For < p < 1, 
the system is degenerate, and many local magnetization 
solutions exist and are reached by hard-spin mean-field 
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FIG. 5: Hysteresis loops for different values of the sweep rate 
h for (a) the pure ferromagnetic phase, p = 0, (b) the fer- 
romagnetic phase with quenched random antiferromagnetic 
bonds, p — 0.15, (c) the spin-glass phase, p — 0.4, all at 
T = 1.5. The loops are, from outer to inner, for sweep rates 
h = 0.05,0.02,0.01,0.005. 

theory. The phase diagram (Fig[T]) is obtained from tem- 
perature T — J -1 and concentration p scans of the equi- 
librium spin-glass order parameter — -j^£jTO? and 
magnetization = -^X^m.;, illustrated in FigJU ob- 

tained by averaging over 20 realizations for aJV = 20 3 
spin system. The results do not change if a larger sys- 
tem is used. In the resulting phase diagram shown in 
Figffl the transition temperatures are gau ged by com- 
paring Tc at p = 0: The precise value [2J] is 4.51, the 
ordinary mean-field value is 6, the value obtained here is 
5.06 . Thus, the transition temperatures are exaggerated 
as expected from a mean-field theory, but considerably 
improved over ordinary mean-field theory. Our obtained 
transition concentrations between the ferromagnetic and 
spin-glass phases are p — 0.22, in excellent agreement 
with the precise value of p — 0.23 pBT ]. 

FiglU shows the zero-temperature spin-glass order pa- 
rameter as a function of antiferromagnetic bond con- 
centration p. It seen that, as soon as frustration is in- 
troduced via the antiferromagnetic bonds, order does not 
saturate at zero temperature, both in the ferromagnetic 
and spin-glass phases, the latter of course showing more 
unsaturation. Moreover, the left column of FigQ] shows 
the equilibrium local magnetizations rrii in a cross-section 
of the system, in the ferromagnetic and spin-glass phases. 
These magnetization cross-sections are remarkably simi- 
lar to the renormalization-group results (26| and are con- 
sistent with the chaotic rescaling picture of the spin-glass 
phase [271 ] . 

Spin- Glass Hysteresis Loops - The quenched random 
magnetic field that is conjugate to the microscopic order 
is Hi{t) — HQ(t)mf^ in Eq.lJTJ, where the are the 
equilibrium local magnetizations obtained with Eq.© 
for a given T,p . Hysteresis loops in the spin-glass or- 
der Q(t) — ■^S,mj(t)m^°' are obtained in the ordered 
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FIG. 6: Spin-glass order parameter Q(t) and uniform magne- 
tization M(t) curves obtained when, in the spin-glass phase, 
the uniform magnetic field H(t) is turned on and then off with 
sweep rate h = 0.005. In this figure, p = 0.4, T = 1.5. 

phases, spin-glass or ferromagnetic, by cycling Hq(t) at 
constant T,p, via a step of magnitude h for each time 
unit. Thus, at time t = 0, Q(t = 0) = Q (0) , the equilib- 
rium spin-glass order parameter. A time unit is N updat- 
ings of Eq.|J5J| at randomly selected sites. Thus, h is the 
sweep rate of the linearly driven [H, 0, Q magnetic field. 
For comparison, it takes of the order of 40jV random up- 
datings to establish equilibrium in the pure ferromagnetic 
phase and 2000iV random updatings to estabish equilib- 
rium in the spin-glass phase. Thus, during the sweeps, 
each time step with N updatings systematically moves 
in the direction of the changed external field but, within 
the loops, does not catch up with equilibrium. The re- 
sulting hysteresis curves are illustrated in FigsO After 
one cycling, the subsequent hysteresis loops for a given 
sweep rate coincide, and are shown in Figs O and used in 
the scaling analysis further below. 

Cycling Effect of a Uniform Magnetic Field on Spin- 
Glass Order - As a contrast to the hysteretic effect of 
the conjugate quenched random magnetic field Hq(t) in- 
troduced above, Fig[6] shows the effect on the spin-glass 
phase of turning on and then off a uniform magnetic field 
H(t) at a sweep rate h. As expected, the spin-glass order 
Q(t) starts at a finite value and returns to zero, while the 
uniform magnetization M(t) — iSjTOj(t) starts at zero 
and returns to a finite value. 

Spin-Glass Hysteresis Area Scaling - The energy dissipa- 
tion of a first-order phase transition is obtained from the 
hysteresis area A of the Q — Hq curve: A — § QdHQ. 
At fixed T,p, the loop area A decreases with decreas- 
ing sweep rate h and finally reaches a value of Ao. The 
area can be scaled as A = A + f(T)h b .^ The (A - A ) 
versus sweep rate h scaling curves are shown in Figs 13 
for the pure ferromagnetic, quenched random-bond ferro- 
magnetic, and spin-glass phases for various temperatures, 
where Aq is fitted. The resulting sweep-rate exponents 
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FIG. 7: The hysteresis area A — Ao versus sweep rate h scal- 
ing curves for T = 1.0(«), 1.5(A), 2.0(4) at different concen- 
trations p. 
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FIG. 8: The sweep-rate scaling exponent b versus concentra- 
tion p for T = 1.0, 1.5, 2.0. These results are obtained by 
averaging over 10 realizations, with the standard deviation 
being used as the error bar. 



b are given in Fig[8]and Table I. From these results, we 
deduce that in the pure ferromagnetic phase, p = 0, the 
exponent b is independent of temperature, as found pre- 
viously Q. However, the value of b = 0.64 that we find 
here, under hard-spin mean-field dynamics, is distinctly 
different from that of b = 2/3 found in Ref.0] under or- 
dinary mean- field dynamics, thereby constituting a dif- 
ferent dynamic universality class. By contrast, in the 
quenched random-bond ferromagnetic phase and in the 
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0.47±0.04 
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0.38±0.01 


0.39±0.02 
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TABLE I: The sweep-rate scaling exponents b at different 
temperatures and concentrations in the ferromagnetic and 
spin-glass phases. 



spin-glass phase, the value of b is distinctly smaller than 
that in the pure ferromagnetic phase, and dependent on 
temperature. Across both of these two phases, there ap- 
pears to be no dependence of b on concentration. 
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